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Several field theoretical approaches to the superconducting phase tran- 
sition are discussed. Emphasis is given to theories of scaling and renor- 
malization group in the context of the Ginzburg-Landau theory and its 
variants. Also discussed is the duality approach, which allows to access 
the strong-coupling limit of the Ginzburg-Landau theory. 

1. Introduction 

The Ginzburg-Landau 
(GL) phenomenological theory of superconductivitjd was proposed long 
before the famous Bardeen-Cooper-Schrieffer (BCS) microscopic theory of 
superconductivity.^ A few years after the BCS theory, Gorkov^ derived 
the GL theory from the BCS theory. This has been done by deriving an 
effective theory for the Cooper pairs valid in the neighborhood of the crit- 
ical temperature T c . The modern derivation proceeds most elegantly via 
functional integrals, introducing a collective quantum field A(x, t) for the 
Cooper pairs via a Hubbard-Stratonovich transformation and integrating 
out the fermions in the partition function.EI 

Amazingly, the GL theory has kept great actuality up to now. It is 
highly relevant for the description of high-T c superconductors, even though 
the original BCS theory is inadequate to treat these materials. The success 
of the GL theory in the study of modern problems of superconductivity 
lies on its universal effective character, the details of the microscopic model 
being unimportant. In the neighborhood of the critical point the GL the- 
ory possesses a wide range of application, many of them outside the field 
of the superconductivity. Perhaps, the most famous example in condensed 
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matter physics is application to the smectic-A to nematic phase transi- 
tion in liquid crystalsl^In elementary particle physics, Gorkov's derivation 
of the GL theory has been imitated starting out from a four-dimensional 
relativistic version of the BCS model, the so-called Nambu-Jona-Lasinio 
model. The GL field in the resulting GL theory describes quark-antiquark 
bound states, the mesons 7r, a, p, and ^iP^This model is still of wide use 
in nuclear physics. Another relativistic GL model is needed to generate the 
masses of the vector bosons W and Z in the unified theory of electromag- 
netic and weak interactions in a renormalizable wayM- 1 This is done by a 
nonabelian analog of the Meissner effect, which in this context is called 
Higgs mechanism. 

In 1973 Coleman and Weinberg^! showed that the four-dimensional 
abelian GL model (or scalar electrodynamics, in the language of elemen- 
tary particle physics) exhibits a spontaneous mass generation in an initially 
massless theory. In the language of statistical mechanics, this implies a first- 
order transition if the mass of the scalar field passes through zero. One year 
later, Halperin, Lubensky, and Ma (HLM) observed a similar phenomenon 
in the three-dimensional GL theory of superconductivity.® These papers 
inaugurated a new era in the study of the GL model. It was the first time 
that renormalization group (RG) methods were used to study the super- 
conducting phase transition. 

At the mean-field level, the GL model exhibits a second-order phase 
transition. The HLM analysis, however, concludes that fluctuations change 
the order of the transition to first. The phase transition would be of second- 
order only if the number of complex components of the order parameter 
is absurdly high. If the number of complex components is N/2, the one- 
loop RG analysis of HLM gave the lower bound N > N c = 365 for a 
second-order transition. No infrared stable charged fixed po int was found 
for N < N c . However, years later Dasgupta and Halperiri-W raised doubts 
on this result by using duality arguments and a Monte Carlo simulation of 
a lattice model in the London limit, to demonstrate that the transition for 
N = 2 in the type II regime is of second order. The RG result that the 
phase transition is always of first order seem s th erefore to be an artifact 
of the e-expansion. Shortly after this, Kleinerll^ perfor med a quantitative 
duality on the lattice and found a disorder field theorjE^ which clarified 
the discrepancy between the HLM RG result and the numerical simulations 
of Dasgupta and Halperin. He showed that there exists a tricritical point 
in the phase diagram of the superconductor on the separation line betw een 
first- and second-order regimes. This result and the numerical predictiorlHI 
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were fully confi rmed only recently by large scale Monte Carlo simulations 
of Mo et a/P3 

In the last ten years the R G approach to the GL model was revisited by 
several groupsff ^ 1 5 | 1 7 | 1 8 | rp-^g a j m was ^ improve the RG analysis in 
such a way as to obtain the charged fixed point for N — 2 predicted by the 
duality scenario. Despite considerable progress, our understanding of this 
problem is far from satisfactory. 

In this paper we shall review the basic field theoretic ideas relevant for 
the understanding of the superconducting phase transition. The reader is 
supposed to have some familiarity with RG theory and duality transfor- 
mations. For a recent review on the RG approach to the GL model which 
focus more on resummation of e-expansion results, including calculations 
of amplitude ratios, see Ref. IT§1 

2. Review of the HLM theory 

The GL Hamiltonian is given by 

n=-(Vx A) 2 + |(V - zeA)^| 2 + m 2 H 2 + ||V| 4 , (1) 

where ip is the order field and A is the fluctuating vector potential. The 
partition function is written in the form of a functional integral as 

Z = J VAV^Vi/j det(- V 2 )<5(V • A) exp (- J d^rT^j , (2) 

where a delta functional in the measure of integration enforces the Coulomb 
gauge V • A = 0. The factor det(— V 2 ) is the associated Faddeev-Popov 
determinant.^ The above Hamiltonian coincides with the euclidian version 
of the Abelian Higgs model in particle physics. The Coulomb gauge is the 
euclidian counterpart of the relativistic Lorentz gauge. The Faddeev-Popov 
determinant should be included in order cancel a contribution 1/ det(— V 2 ) 
that arises upon integration over A taking into account the constraint V • 

a = o.EEOl 

2.1. HLM mean-field theory 

Let us make the following change of variables in the partition function J5J : 
yj = l=pe* e , A = a+ivf?, (3) 
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which brings the partition function to the form 



Z = JvaVpVOp det(- V 2 )<5(V • a + e" 1 V 2 6») cxp (- J d 3 rH^j , (4) 
with 

H = \{V x a) 2 + ^p 2 a 2 + \{V P f + ^p 2 + (5) 

In these variables, the Hamiltonian does not depend on 0. This change of 
variables is allowed only in a region where the system has few vortex lines, 
since otherwise the cyclic nature of the ^-ficld becomes relevant, and V e 10 
is no longer equal to iVO e 10 , as assumed in going from Q to JSJ, but equal 
to z(V0 — 27rn)e ie , where n is a vortex gauge field. I -^ I ^ I This problem 
of the Hamiltonian JSJ which is said to be in the unitary gauge will need 
special attention. 

Now we can use the delta function to integrate out 9. The result of this 
integration cancels out the Faddeev-Popov determinant. 

Next we assume that p is uniform, say p = p = const. Since the Hamil- 
tonian is quadratic in a, the a integration can be done straightforwardly to 
obtain the free energy density: 

^=^Trln[(-V 2 + e 2 p 2 )^ + ^]-^^(0)Trln( e V) + ^p 2 + ^- 4 , 

(6) 

where V is the (infinite) volume and the term J 3 (0)Tr ln(e 2 p 2 )/2y comes 
from the exponentiation of the functional Jacobian in Eq. Q. In order to 
evaluate the Trln in Eq. we have to use the Fourier transform of the 
operator 



M(r, r') = [(-V 2 + e 2 p 2 )S^ + d^}6 3 (r - r'). (7) 
The operator M(r, r') is diagonal in momentum space. Indeed, we can write 

M(r,r') = J ^ e *P-(-r')M(p), (8) 

M(p) = (p 2 + eV)i^(p) + e 2 p 2 P^(p), (9) 

where P^ v {p) = <^„ - P^Pu/p 2 and P^{p) = P^Pv/p 2 are the transverse 
and longitudinal projectors, respectively. Now, the Trln M(r, r') is obtained 
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first by taking the logarithm of the transversal and longitudinal parts of 
A/(p) and tracing over the vector indices. The second step is to trace out 
the coordinates by integrating over r for r = r'. This produces an overall 
volume factor V which cancels out exactly the 1/V factor in Eq. JJjJl. Note 
that the term proportional to S 3 (Q) ln(e 2 p 2 ) drops out in the calculations. 
If we write p = \tf>\, the end result is the celebrated HLM mean field free 
energy: 

^H L M = -^ + ^M 2 + ||#. (10) 

Due to the cubic term in Eq. (|1U|I . the transition is of first order. The basic 
problem with this argument was pointed out in Ref. 11 In the critical 
regime of a type-II superconductor, the order field contains numerous lines 
of zeros which make it impossible to use uniformity assumption p — p = 
const. 

2.2. Renormalization group in d — 4 — e dimensions 

In the original HLM work, the RG calculations were performed using the 
Wilson version of the renormalization groupp^ where fast modes are inte- 
grated out to obtain an effective theory in terms of the slow modes. The field 
theoretical approach using the Callan-Symanzik equation gives an equiva- 
lent result in a perturbative setting. We shall discuss the RG calculation in 
d — 4 — e dimensions in this context. 

The dimensionless renormalized couplings are defined by 

f = p-'Z A el, g=u-'ZlZ- x u , (11) 

where p gives the mass scale of the problem. Note that in Eq. (|llf> we have 
denoted the bare couplings by a zero subindex. The bare fields are denoted 
by -00 and Ao. Accordingly, the bare mass is denoted by mo- We shall use 
this notation from now on. 

The renormalization constants are defined such that the "renormalized 
Hamiltonian" Ti. r is given by the following rewriting of the bare Hamilto- 
nian: 

H r (A,V;m 2 , W ,e) =H{Z 1 fA,Z 1 J 2 ^-Z m Z- 1 m\Z g Z- 2 u,Z^ 1/2 e). (12) 

From Eq. Ijl2|l we see that the renormalized fields are given by ip — Z,^ 1 ipo 
and A = Z A ^ 2 A . 
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The calculations are more easily done if we set m — and evaluate 
the Feynman graphs at nonzero external momenta to avoid infrared di- 
vergences. The external momenta in a four leg graph will be taken at the 
symmetry point: 



Pi- Pj = ^(4*y-1). (13) 

Even if to ^ we have to face severe infrared divergences in this prob- 
lem. This is due to the masslessness of the vector potential field. The free 
A propagator is given in the Coulomb gauge by 



^(P) = ^(V-^). (14) 
Thus, the graph in Figgis proportional to the integral 

d d k 1 



(27r) d k 2 (p-k) 2 



tt^ 2 r(2-d/2)r 2 (d/2-l]^ 



SP 



(2?r) d r(d-2) 



evaluated at the symmetry point (SP) as prescribed in Eq. 11311 . 





Fig. 1. Feynman graph contributing to the / term in fi g . The dashed lines represent 
the vector potential propagator, while the external lines represent the order parameter 



The /3-functions are defined by 

0/ = ^. & = >v (16) 

These /3-functions are defined for the dimension d in the interval (2, 4]. For 
d > 4 the theory is no longer renormalizable and the /3-functions are not 
defined. At d = 2 the infrared divergences are very strong and require a 
special treatment. 
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Let us assume that the order parameter field has N/2 complex com- 
ponents. Then, for any dimension d G (2,4] the /3-functions are given at 
one-loop order by 

f3 f = (4-d){-f + NA(d)f}, (17) 



P g = {A~d){ -g + B(d) 
where 



-2(d-l)/.9+^±V + 2(d-l)/ 2 



(18) 



If we set d = 4 — e and expand to first order in e, we obtairPH 

^/ = -/ +i ^/ 2 , (2D 

From Eqs. (|21|l and (12211 we see easily that charged fixed points exist only 
if N > N c — 365.9. Thus, if N is large enough to allow for the exis- 
tence of charged fixed points we obtain the flow diagram shown schemat- 
ically in Fig. |3 In the figure the arrows correspond to /x — ► 0. There are 
four fixed points. The Gaussian fixed point is trivial and corresponds to 
f+ = g„ =0. It governs the ordinary mean-field behavior. There is one 
non-trivial uncharged fixed point, labeled "Heisenberg" in the figure, which 
governs the iV-component Heisenberg model universality class. For N = 2 
the superfluid 4 He belongs to this universality class (in this case we speak 
of a XY universality class) . The Heisenberg fixed point is unstable for non- 
zero charge. There are two charged fixed points. The one labeled SC in the 
figure is infrared stable and governs the superconducting phase transition. 
Its infrared stability ensures that the phase transition is second-order. The 
second charged fixed point is labeled with a T and is called the tricritical 
fixed point. It is infrared stable along the line starting in the Gaussian fixed 
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point and unstable along the g-direction. The line of stability of the tricrit- 
ical fixed point is called the tricritical line. The tricritical line separates the 
regions in the flow diagram corresponding to first- and second-order phase 
transition. 



As we shall see later on in these lectures, the flow diagram shown in 
Fig. |21 should also be valid for N = 2. a For the moment let us remark 
that if we use the /^-functions in fixed dimension as given in Eqs. Ijl7(l and 
(I18fl the value of N c is considerably smaller at d = 3. Indeed, by setting 
d = 3 in Eqs. (|17(l and . we obtain that charged fixed points exist for 
N > N c = 103.4. 

2.3. Critical exponents 

The critical exponents can be evaluated using standard methodsP2 The i] 
exponent is obtained as the value of the RG function 



f 




Heisenberg 



Fig. 2. Schematic flow diagram for the GL model for N > N c . 



7V> = M 



(23) 



a In the context of the e-cxDansion, sec for example the two-loops resummed RG analysis 
by Folk and Holovatch. Ed 
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at the infrared stable fixed point. The rj exponent governs the large distance 
behavior of the order field correlation function at the critical point: 

(#^t (r0) ^ (24) 

The critical exponent v, governing the scaling of the correlation length 
£ = mT 1 ~ t~ u , with t being the reduced temperature, is obtained as 
the infrared stable fixed point value of the RG function defined by the 
equation 



— - 2 = 7 m, (25) 

v-4, 

where 



^=4Xi)- (26) 

In an approach where the correlation functions are computed at the critical 
point the mass renormalization Z m must be computed through |?/>| 2 inser- 
tions in the 2-point functionP2E For any dimension d 6 (2,4], we obtain 
the one-loop result: 

^ = (l-d)(A-d)B(d)f, (27) 



7m = (N + 2)(d - i)B(d)g/2 - (28) 

Once the critical exponents rj and v are evaluated, all the other critical 
exponents can be obtained using the standard scaling relations. We shall 
prove later that the standard scaling relations apply also in the case of the 
GL model. 



2.4. 1/N expansion 

The 1/N expansion is one of the most popular non-perturbative methods in 
the field theoretical and statistical physics literature. The critical exponents 
of the GL model at 0(1/N) were calculated at d = 3 by HLM. Since the 
original paper does not contain any detail of the calculation, we shall outline 
it here. The large N limit in the GL model is taken at Nu and Ne 2 fixed. 
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Let us calculate the critical exponent r\. The relevant graphs are shown in 
Fig-El The idea is to pick up the p 2 ln|p| contribution of r^ 2 )(p) = G _1 (p), 
where G is the order field propagator. To keep the theory massless, we 
have to subtract the p = contribution of the self-energy. The subtracted 
contribution coming from the graph (a) in Fig. is 



S«(p) 



d 3 k 



k 2 



where 



V(k) = 



2 a 



1 + «JVn(k) ' 



and n(k) is the polarization bubble: 

, d 3 k 

n(p) 



G (p-k)G (k), 



(29) 



(30) 



(31) 



where Gq is the free scalar propagator. In the massless case we have simply 



n( P ) 



i 

8|p| 



(32) 



When e = Eq. (|30(l gives the effective interaction of the O(N) modePl 
The —rjaP 2 In |p| contribution from Eq. I|29[l gives 



(33) 



3tt 2 ^' 

which is just the ?y-exponent of the O(N) model at 0(l/./V)Pl 

In order to compute the contribution from the graph (b) of Fig. 10, 
we need to calculate the dressed vector potential propagator. The dressed 
propagator is obtained by summing the chain of bubbles shown in Fig. [31 
The result is 



£V(p) = 



l 



P 2 + Ea(p) 1 

where in the massless case 



5n U - 



(34) 



Ne 2 

£4p) = -^-IpI 



(35) 
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Ca) 




Fig. 3. Feynman graph contributing to the momentum dependent part of the self-energy 
of the order field propagator. The double lines represent dressed propagators. The graphs 
(a) and (b) contain effective vertices given by an infinite sum of a chain of bubbles. 



The self-energy contribution from graph (b) is then 



S b (p) 



d 3 k (2p IM -k IM ){2p u -k u ) 



(2tt) 3 (p-k) 2 
Thus, we obtain the contribution — r/bp 2 In |p|, where 



(36) 



r]b = 



128 

'3n 2 N 



(37) 



The 77-exponent of the GL model at 0(1/N) is obtained by adding the 
contributions from Eqs. (|33|l and l|37|l : 



W = Va + Vb 



40 
n 2 N' 



(38) 



The critical exponent v at 0(1/ N) can be computed through similar 
calculations, except that one needs to consider massive propagators for 
the order field. Instead picking up the contribution p 2 ln|p|, we take the 
m 2 In m one from the self-energy at zero momentum. This gives in fact the 
critical exponent 7. The critical exponent v is obtained straightforwardly 
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by using the scaling relation 7 = v[2 — 77). The end result is: 



-'-as?- < 39 > 

Note that in the 1 /N expansion no tricritical fixed point is generated at 
0(1 /N). The reason for that comes from the fact that the graph in Fig. ^ 
is 0(1 /TV 2 ) and does not contribute to the above calculations. This graph 
is the main obstacle against attaining the charged fixed point. 



3. Existence of the charged fixed point 
3.1. Scaling near the charged fixed point 

By assuming the existence of the charged fixed point, it is possible to derive 
exact scaling relations for the GL theory. 1 ^ 21 For example, it is easy to prove 
that if a charged fixed point exists then v' — v 1 where v' is the exponent 
of the penetration depth A. To this end, it is necessary to consider the GL 
model in the ordered phase, that is, at T < T c . In this situation the vec- 
tor potential becomes massive through the Higgs mechanism.pl with a mass 
rriA = A -1 . The Higgs mechanism implies that there are no massless modes 
in the superconductor in the ordered phase. This shows a fundamental dif- 
ference between a superconductor and a superfluid: The ordered phase of 
a superfluid contains a massless mode or Goldstone boson. The Ginzburg 
parameter k is defined by the ratio between the Higgs mass and the vec- 
tor potential mass: k = m/mA- It can be shown that m? — up s /2 and 
m 2 A = e 2 p s , with p s being the superfluid density. These formulas are easily 
obtained in a tree level analysis of the Higgs mechanism. Their application 
in the renormalized case follows f rom im posing suitable renormalization 
conditions and the Ward identities.^SOl Thus, k 2 = u/2e 2 = g/2f. Since 
both / and g tend to a nonzero fixed point value in the infrared (m — * 0), 
it follows that the Ginzburg parameter k 2 — A 2 /£ 2 = g/2f — > const as 
m — > 0. As a consequence, the scaling relation A ~ £ holds and therefore 
the equality between the corresponding critical exponents. 

In the following it will be useful to use m as the running RG scale. Thus, 
/ = m d - i Z A el and g = m d - 4 Z 2 Z g - 1 u . Therefore, 



df 

P f = m^ = ( 7A +d-A)f, (40) 
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where 

dlnZ A 

1a = m gm . (41) 

Thus, under the assumption of existence of the charged fixed point we 
obtain that 

r) A =7A(f*,9*)=4-d. (42) 

The abo ve equ ation gives the exact anomalous dimension of the vector 
potential! I ' This means that the large distance behavior at the critical 
point: 

(Mr)Mr')) ~ |r _ r1 L + , A - j^^Tp (43) 

which holds for all d € (2,4]. 

It is instructive and experimentally relevant to compare the scaling near 
a charged fixed point with the one governed by the XY fixed point. To this 
end, let us first note that the flow equation for n 2 is written exactly a 

m ^L = K * + 4 - d - 7 A (44) 
om \ g J 

Since k 2 — m 2 /m\, it follows the following exact evolution equation for 



m^ = m^d-2 + lA -^y (45) 
From Eq. I|45|l we obtain that near a fixed point via behaves as 



m A ~ m ( d -i+VA)/^ ( 46 ) 

The abo ve c onstitutes a rederivation of a result due to Herbut and 
Tesanovic^l For the charged fixed point t]a — 4 — d and we obtain once 
more that v' = v (remember: m = £ _1 and tua — A ~ 1 ). For the XY fixed 
point, on the other hand, t\a = and we obtain from Eq. Q46[) the scaling 
relation for the superconducting XY universality classPSl 

,' = ^. (47) 
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A further interesting consequence of the scaling relation (|46|l is the fol- 
lowing. The vector potential mass is given in terms of the superfluid density 
p s as m 2 A = e 2 p s . Eq. H4U|) implies the scaling e 2 ~ m VA near the charged 
fixed point. Therefore, from Eq. I)46|l we obtairP^ 

p s ~ m d - 2 ~ \t\ v ( d - 2 \ (48) 

where t is the reduced temperature. Eq. (|48|l is just the Josephson relational 
and we have thus shown that it is also valid near the charged fixed point. 
In his original paper, Josephson has obtained the above relation for the 
superfluid in the form p s ~ t 2 0~ v v^M Then he derived the result (|4^|) by 
assuming that hypcrscaling holds, that is, dv = 2 — a, which together with 
the scaling relations a+2/3+7 = 2 and 7 = 2/(2— ?y) imply 2(3— vri = v(d—2). 
We have obtained the result without using these supplementary scaling 
relations. Our result follows from the exact evolution equation for the vector 
potential mass, Eq. g5j). The form p s ~ t 2f) ~^ can be proven without any 
reference to the charge and is therefore valid also for the superconductor. 
From this statement and Eq. (|48J) we prove 

2/3 - vq = v(d - 2) (49) 

for the superconductor. 

Experiments in high quality crystals of YBa2CusOj-s (YB CO) per- 
formed at zero field verify very well Eq. I|47(l for the d — 3 caseM^ Most 
experiments are unable to probe the charged critical region. However, in ex- 
periments involving critical dynamics the situation is not clear. In this case 
we have again that different scaling relations are obtained near the charged 
fixed point. In general the AC conductivity scales as cr(w) ~ e 2 p s /(— ioj). 
Since p s ~ £ 2 ~ d , e 2 ~ £,~ VA , and u> ~ w here z is the dynamical critical 
exponent, we derive the scaling relation^ 

For the XY universality class where T]a — we have o~(u>) ~ £ 2 ~ d+z .E 3 
Near the charged fixed point, however, we obtain cr(uj) ~ |^|^(2-2)) 28 | 29 | 

3.2. Duality 

Duality is a powerful tool in physics. ^ It allows the mapping of a weak 
coupling problem on a strong-coupling one. In the context of statistical 
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physics, it maps the low temperature expansion into a high temperature 
expansion. In some cases, duality allows to obtain exact information on the 
physical system. The classical example is the two- dim ensional Ising model, 
where the exact critical temperature was obtained^ before the exact solu- 
tion appeared^ The exact determination of the critical temperature was 
possible because the Ising model has the self-duality property, that is, the 
duality transformation has a fixed point. The self-duality property is also 
verified in other systems and for d > 3, like in the Zi lattice gauge the- 
ories. However, the discreteness of the gauge group makes these theories 
very similar from the point of view of duality to the two dimensional Ising 
model. Self-duality is more difficult to be found in continuous gauge groups. 
The GL model, for example, has no such property, but as we shall see, it is 
nevertheless almost self-dual. 

In this section we shall discuss the field theoretical approach to duality 
in the GL model. We shall show in detail how scaling works in a disorder 
field theory (DFT) for the superconducting phase transition. The DFT 
to be discussed here was proposed first by Kleinert nearly twenty years 
a go_El This formulation allowed to demonstrate that a tricritical point 
exists in the phase diagram of the superconductor. The existence of this 
tricritical point allowed to build a consistent picture where the strong- 
coupling limit - which exhibits a second-order phase transitional - and 
the weak coupling limit, with its weak first order scenario, coexist with 
the normal phase, meeting at the tricritical point. On the basis of the 
DFT, the estimated value of k at the tricritical point was n t ~ 
Early Monte Carlo simulationJ^Hl gives, on the other hand, the estimate 
Kt ~ 0.42/v2- Remarkably, Kleinert's estimate agree within 5 % with a 
recent, more precise, Monte Carlo simulation by Mo, Hove, and 

Sudb0PI 

The first scaling analysis of the DFT was made by Kiometzis, Kleinert, 
and SchakelPI From the analysis in Ref.1141 it was possible to establish the 
value of the critical exponent v as having a XY value, v ~ 0.67. However, 
as we shall see, the scaling analysis of the DFT contains some ambiguities 
which are not yet completely solved. 

3.2.1. Duality in the lattice Ginzburg- Landau model 
A lattice version of the GL model has the Hamiltonian 



H = -/3^ cos(V^ ~ eA ^) + 2 E( V x A *) 2 ' 



(51) 
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where is the lattice derivative, V^/j = — /.;, and (3 = 1/T. The 
partition function is then given by 



Z 



n 



2tt 



i,fj 



exp(-F). 



(52) 



The duality transformation can be done exactly w hen the Villain form of 
the Hamiltonian is used. The Villain approximation^^ corresponds to the 
replacement 



o cos x 



E 



e 2 



(x-2-rm) 



(53) 



n— — oo 



which turns out to be very accurate near the critical region!*^ 

Assuming from now on the Villain approximation, we introduce an aux- 
iliary integer field such that 



E exp 







oc exp 



(54) 



The proportionality factor above is not important in the following. All such 
proportionality factors will be neglected in the foregoing manipulations. 
They correspond to smooth factors in the temperature. Eq. 1)54(1 was ob- 
tained using the identity 



E eM 



/2)m -\-ixm 



= \ - E e(_1 



/2t)(x-27rn) 2 



(55) 



The summation notation in Eq. 1|54J) with a {rii^} means a multiple sum- 
mation, analogous to multiple integration. 

By integrating out the angular variables 9i we obtain the partition func- 
tion: 
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Z 



E ^V mi.oexp I ^2 

{m,} I i 



2/3 



m~ + ieAi ■ m, 



-(V x Aj) 2 



(56) 



The Kronecker delta constraint V • ir^ = generated by the Oi integrations 
implies that the links variables form closed loops. These closed loops 
arc interpreted as magnetic vorticesP^ When e = the vector potential 
decouples and we have, up to proportionality factor, the partition function 
for the XY model in terms of link variables: 



Zxy = <5 v . mi ,oexp 



(57) 



We can solve the constraint on mi by introducing a new integer link 
variable through = V x Ij. After integrating out A^ and using the 
Poisson formula 



F(n) = / dxF(x)e 27rmx (58) 

n— — oo m— — oo — 00 

to go from the integer variables h to continuum variables hi, we obtain 



Z 



Y[dhi 



fm,} 



E 



2/3 



(V x h, 



— h? + 27rim. J • hi 



(59) 



Eq. H59|) corresponds to the dually transformed lattice GL model. 
By performing the integration in Eq. I|59|) we obtain 



Z = 2J <5v- mi ,oexp 



-2tt 2 /3 ^ m ia G{Yi - Yj)m jtl , (60) 
where the Green function G has the following behavior at large distances: 

g-V/Jelri-rjl 



G(r, - r,) 



47r|r 4 



(61) 
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Thus, in the superconductor the magnetic vortex loops interact with a 
screened long range interaction. 

Let us consider now the "frozen" superconductor limilP^ of the dual 
representation (|59|) . The "frozen" superconductor corresponds to the zero 
temperature limit, T — > 0. In this case, after integrating out h^, we obtain 



-^frozen — ^V m i: ex P ( jj - 




(62) 



Eq. (|62|l has the same form as Eq. I|57|) . Thus, the "frozen" superconductor 
is the same as a XY model provided we identify 



4tt2 
T 



(63) 



The above result allows us to localize a critical point on the e 2 -axis in the 
phase diagram in the e 2 — T-plane. Using Eq. (|63l) and the fact that T c «3 
for the XY model in the Villain approximationj^we obtain e 2 = An 2 /T c k, 
13.159 on the e 2 -axis. Thus, we can locate two limiting critical points in the 
phase diagram, since we have in the T-axis the Villain-XY" critical point at 
T c «3. 

The vortex-vortex interaction in Eq. H6Q(I is singular at short distances. 
Therefore, it is natural to introduce a vortex core term with energy eo m 
the dual lattice Hamiltonian: 



E 



^(V x h,) 2 + yh 2 - 27u mi ■ h, + |m 2 



(64) 



Using Eq. (|55|) and the integral representation 
of the Kroncckcr delta 



we obtain 



JVmi,0 — 



2ir 



(65) 



— (V x h,) 2 + -h 2 

2/3 v ' 2 



2itin 



2tt^) 5 



(66) 
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The Hamiltonian in Eq. I|6(jf> has the same form as the lattice GL Hamil- 
tonian in the Villain approximation, except that in Eq. I|66[l the vector 
field is massive. Thus, we see that there is almost a self-duality between 
them. When e = the Hamiltonian (|66|) is the dual of the (Villain) lattice 
XY Hamiltonian. Note that in this duality transformation a locally gauge 
invariant model is mapped on a globally invariant one. 

Let us set e = in Eq. Ij64|l and look for the phase diagram in the j3 — eo 
plane. In such a phase diagram the point (/3 C , 0) corresponds to the XY 
critical point. Integrating out we obtain the partition function: 



Z'\ e=0 = X <5v. m< ,oexp 



-2tt 2 /3 X m iA1 G(ri - r,)m iM - y m 2 



, (67) 



where G behaves at large distances as 



From Eqs. I|57|) and Ij67|l we see that the point (0, l/20r) in the (3 — eo plane 
corresponds to an "inverted" XY (IXY) transition!"^" By performing the 
Ai integration in Eq. (|56|) we obtain 



Z 



E 

in..} 



<Vm,,oexp 



2/3 



(69) 



and we see that the IXY critical point corresponds to e 2 = 47r 2 /3 c .E3 

Clearly the IXY transition is a second-order phase transition. Note that 
this transition arises in a lattice GL model where the amplitude fluctuations 
are frozen (London limit). Therefore, we should not expect to find a first- 
order phase transition in this case. This London limit is appropriate when 
magnetic fluctuations are relevant in the type II regime. In Fig. 0] we the 
approximate phase diagram.UJ 



3.2.2. The disorder field theory 

The IXY universality class must have the same thermodynamic exponents 
as the XY model. For instance, we expect v ~ 0.67 and, from the scaling 
analysis of Section 3.1, v 1 = v. This last scaling relation has been confirmed 
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Fig. 4. 
Here T c 



Schematic phase diagram showing the critical points on the e 2 - and T-axis. 
s; 3. The ordered superconducting phase corresponds to < e 2 < e 2 . 



in Monte Carlo simulations of the lattice model i|51|) .^"3 Although the ther- 
modynamic exponents are the same as in the XY universality class, critical 
exponents as rj and t\a are not the same. We have seen in Section 3.1 that 
tja = 4 — d near the charged fixed point. The d = 3 value, t]a — 1, corre- 
sponds to the IXY universality class discussed in the preceding subsection. 
The XY universality class, on the other hand, has tja = 0. Also, we have 
— 1 < rj < in the IXY universality class, while rj > in the XY one. 

Useful information from different universality classes and crossovers in 
superconductors can be obtained from the DFT. The DFT is constructed 
out the dual lattice Hamiltonians discussed in this s ection. The bare DFT 
associated to the lattice Hamiltonian 1)66(1 is given b y^ l ^ l 



Hdft 



J(Vxh ) 2 - 



l 



(V-ieoho^'+m^op+^ol 4 , (70) 



where the bare dual charge eg = 2irmAfi/ e o- ^A.o is the bare mass of the 
vector potential in the original theory. We see that the disordered phase 
of the DFT corresponds to the ordered phase of the GL model. This is a 
general feature of all duality transformations: the low temperature phase 
is mapped in the high temperature phase. This justifies the denomination 
"disorder field theory" for the continuum limit of the lattice dual model. The 
field (f>o is the bare disorder parameter field. In the superconducting phase 
(</>o) = 0, while the order parameter in the original GL model (■0o) 7^ 0. 
Conversely, the normal phase corresponds to (^o) ^ and (tpo) = 0. 

It should be noted that the Hamiltonian (|7t)|l is a generalization of the 
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London model. The field ho is in fact the magnetic induction, while |<^o| 2 
gives the vortex loop density. 

The renormalization of the Hamiltonian l|7U|) is similar to the one of 
the GL model, up to the following subtlety. From the Ward identities we 
obtain that the mass term for the vector field is not renormalized, that is, 
m 2 A h 2 /2 = 77i^h 2 /2, where the absence of the zeroes subindices indicate 
renormalized quantities. Since the renormalized induction field is given by 

— 1/2 

h = Z h ' ho, we obtain 

m\ = Z h m 2 A0 . (71) 
The dual charge renormalizes in a similar way: 

e 2 = Z h e 2 . (72) 

Since e 2 = 47r 2 m 2 4 / e o> ^ follows from Eqs. (fTTJl and fT^ that the Cooper 
pair charge eo is not renormalized in the DFT, e = eo. 

It is important to remark that the vector potential mass renormalization 
in the DFT involves only one renormalization constant, while the same is 
not true in the GL model.^1 

Due to the presence of a massive vector field, the DFT h as a ambig uous 
scaling!^ This ambiguity was a matter of debate recently EZEEESTLet'g 
see how it works. In Ref. I14l the scaling chosen was in principle very natural, 
with the bare masses behaving in the same way with mean field exponents: 
m 2 , ~ \t\ and m 2 A ~ \t\. Renormalization is employed as usual. We have, 

ml = Zf^Z^rh 2 , (73) 

and therefore, 



™^ = (2 + 7m)^o> (74) 
dm 

where 7™ is defined in a way similar to j m in Eq. (|28|) . Since m\jm A = 
const, we obtain 



fh = (2 + 7m) m i,o- ( 75 ) 
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Let us define the dimensionless renormalized dual coupling by / = e 2 /m. 
From Eqs. (|72|l and l|75(l . we obtain the /3- function: 



- df 
Pf = " 

where 



/?, = m-^r = {lh + 7m + 1)/, (76) 
1 am 



_ d\iaZ h 

lh = m ^ m - (77) 

Now, we can easily obtain the following bound 

7A + !<--!< 1. (78) 
v 

It is also straightforward to show that 7/, > 0. Therefore, the infrared stable 
fixed point to l|76|l is at /* = 0. This implies that the critical exponent v 
has a XY value and that r)h = 7^ = 0. 
From Eqs. l(7T|) and J75J we obtain 



m 



dm 

Near the fixed point the above equation becomes 



(7h+7m + 2)mi. (79) 



which implies the scaling 



fn^ « -mi, (80) 



mi - m 1 ^. (81) 



Since to ~ \t\ v the above scaling implies that the penetration depth expo- 
nent is given exactly by v' = 1/2. 

Another possible scaling is the one considered by HerbulP^I where it is 
assumed that m\ = const. Within this scaling we obtain instead Eq. (|76() 
the following /3-function: 



(3f = (lh - 1)1 (82) 
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which is similar to the /3-function of the coupling / in the GL model at 
d = 3. From Eq. f?TJ) we obtain 



= 7fc"4 (83) 

The /3-function for the coupling g = - u /fh contains functions of the ratio 
m/mA multiplying every power of J jsiZI Due to the evolution equation JHUJl, 
we see that m A ~ m ih ~ rh. Therefore, rh/iriA — > as the critical point is 
approached. Thus, fixed point g* is the same a s in . the XY model and once 
more the critical exponent v has a XY value. "*3 However, from the exact 
scaling behavior m 2 A ~ m we see that the penetration depth exponent 
is given by v' — v/2, which corresponds to the same value as in the 3D 
XY superconducting universality class. Therefore, the scaling considered 
in Ref. [22 does not give the expected value for the IXY universality class, 
which should be v 1 = v w 2/3.121 It was shown in Ref. E3 that the IXY 
universality class can only be obtained if the bare mass m 2 A ~ |tp, where 
C = 2v « 4/3, i.e., Tn A Q should scale as m of the original GL model. It 
seems that if we want that the dual model implies the result v' = v, we 
have to make this assumption. 

An alternative scenario for the scaling behavior in the DFT is the fol- 
lowing. We have shown that the XY model dualizes in a GL model. The 
GL model, on the other hand, dualizes on the DFT whose Hamiltonian is 
given in Eq. I|7U|) . Now, the dual of the dual must be of course the original 
model. This means that the GL model should dualize in a XY model and 
therefore the DFT Hamiltonian should be equivalent to it. On the basis of 
this argument we should expect a scaling consistent with the XY univer- 
sality class, instead of IXY . If we accept this argument, we are led to the 
conclusion that the correct scaling behavior should assume m 2 A = const 
as in Ref. 1371 to obtain the XY scaling of the penetration depth, v' = v/2. 



4. The physical meaning of the critical exponent rj 

In the superconducting phase transition only the exponents v, is', and a 
are measured. Here a is the specific heat exponent, which is related to v 
by the hyperscaling relation dv = 2 — a. At present the critical exponent 
X] is not measured and we can even doubt of its physical significance. We 
can argue that the superconducting order parameter cannot be considered 
to be a physical measurable quantity because its conjugate field has no 
physical meaning. In a ferromagnet the field conjugate to the magnetization 
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is just the external magnetic field, which can be controlled by experiments. 
Another problem is that the order parameter (ip) is not gauge invariant. 
Thus, a calculation of (ip(r)ip> (r')) will depend on the gauge choice and, as 
a consequence, r\ will also be gauge dependent. 

In this Section we shall show that it is possible to define a gauge- 
independent rj exponent and discuss its physical significance. The physical 
meaning of r\ arises due to a special feature which at first glance looks very 
much like a pathology: it has a negative sign. Indeed, we would expect from 
very general non-perturbative arguments that 77 should be positive. In the 
case of a pure |i/>| 4 theory we can prove that r\ > in the following way. In 
momentum space the correlation function G(r — r') = (^(r)^ (r')) has the 
spectral representationE^I 

^ (P)= f^^' (84) 
where the spectral weight p(p) satisfies the sum rule 



dfip(n) = 1. (85) 

The above representation is well known in quantum field theory and is called 
Kallen-Lehmann spectral representation.ESI Now, because of the condition 
(|85fl on the spectral weight, we have the inequality 

G(p) < -L. (86) 

If one assumes the low momentum behavior at the critical point G(p) ~ 
l/|p| 2 ~'', we obtain from the inequality l|8fi[) that 77 > 0. 

In the case of the GL model, all c alculations of 77 give a negative value 
in the interval — 1 < i] < for d = 3. ) l 1 $ l 1 v | r t ql 4 vrMl 1 s? l l n general it is 
argued that since G(p) ~ 1 / 1 p> | 2 77 , we have in real space the large distance 
behavior at the critical point: 



G ( r -0~ |r _ r ^ 2+tr (87) 

The above will not diverge as |r — r'| — > 00 provided 77 > 2 — d and for 
this reason we could have in principle a negative rj exponent. Such an ar- 
gument is certainly not correct in the case of pure |^;| 4 theory where the 
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Kallcn-Lehmann representation holds which implies rj > 0. In the case of 
the GL model the situation is much more subtle and the Ka llen-Le hmann 
representation does not apply, at least not in the above formPHBSI 

In order to give a physical interpretation to the negative sign of T] in 
superconductors, let us consider a one-loop approximation at the critical 
point and fixed dimensionality d — 3. The calculation is uncontrolled but 
serves to illustrate the main point. Assuming N = 2, the vector potential 
propagator is then given by 



n ( \ - 1 f A _ PuPv 

p 2 + e 2 |p|/16 \ 

while the order parameter two-point correlation function is 

G (P) = 2 \ 1/7 • (89) 
\> l — e^|p|/4 

Now, we see from Eq. (|88|) that as |p| — > the second term in the de- 
nominator dominates implying t)a = 1. However, the same argument does 
not apply to Eq. I|89|l because the second term in the denominator has a 
negative sign in front of it and therefore the p 2 term is still relevant. Thus, 
there is a momentum space instability in the problem similar to the one 
encountered in theories of magnetic systems exhibiting a Lifshitz point E01 
There the Hamiltonian already contains the momentum space instability 
from the very beginning due to the presence of higher order derivatives. 
E01 D ue to this, the susceptibility in those magnetic systems has a maxi- 
mum at a nonzero value of p. This lead to the appearance of a modulated 
regime in the phase diagram, which is plotted in the P — T-plane, where 
P is the ratio between two competing interactions. It was conjectured in 
Ref. 1281 that a similar behavior occurs in the superconductor. The modu- 
lated regime would be associated to the type II behavior, which would be 
in this way analogous to the helical phase in magnetic systems exhibiting 
a Lifshitz point. In the case of magnetic systems, the Lifshitz point is the 
point in the phase diagram where the paramagnetic, ferromagnetic, and 
helical phases coexist. In the case of the superconductor it corresponds to 
the point where the type I, type II, and normal phases coexist. In straight 
analogy with magnetic systems, the phase diagram is plotted in the k 2 — T- 
plane. Note that in the case of the GL model the Lifshitz point-like behavior 
would be generated by thermal fluctuations. 

Further insight in this problem can be obtained by looking at the propa- 
gator in the 1/N expansion already discussed in this review. The self-energy 
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at the critical point and 0(1/ N) is given by 



Thus, besides the pole at p 



S(p) 



we have also a pole at 




(90) 




7T 



■ 2 N 



) 



(91) 



40 



This instability is similar to the on e le ading to chiral symmetry breaking 
in three-dimensional QED (QED3).EUThe difference is that in QED3 the 
instability occurs with respect to the mass, which is dynamically generated 
by spontaneous chiral symmetry breaking. Thus , if M is the generated 
fermion mass, the pole of the propagator occurs aJ^ 



We have not yet discussed how to cure the disease which results from 
the lack of gauge invariance of the order parameter correlation function 
G(r — r'). Our physical interpretation of the critical exponent r] given above 
has no value if 77 is a gauge-dependent quantity. The best thing to do is to 
define a gauge-invariant correlation function for the order parameter. The 
choice of such a gauge-invariant correlation function is not unique but as 
we shall show, there is one whose value of r\ coincides with the one which 
is obtained by computing it in the Coulomb gauge, which is the gauge we 
are using in this review. 

A popular gauge-invariant correlation function which is often used in 
the literature is 



A calculation of rj using the ab ove correlation function was carried out 
recently by Kleinert and Schakel!^ They calculated this exponent using 
both the e-expansion and the 1/A-expansion. In the former case the result 
is 




(92) 




(93) 



36 



(94) 
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while in the latter case 77 was computed for arbitrary dimensionality d £ 
(2,4) and up to order 1/A^2 

__ 1 (d 2 + 2d-6)T(d-2) 



' NT(2~d/2)T 2 (d/2- l)r(d/2)" 

By setting d = 4 — e in Eq. (|95|l and expanding to order e we obtain Eq. 
(|94|l . In a GL Hamiltonian with a gauge- fixing term 

^gf-T^(V-A) 2 , (96) 

the above expressions of 77 corresponds to values that would have been 
obtained by fixing the gauge a = — 3 in the e-expansion case and a = 1 — d 
in the case of the l/TV-expansion. Thus, in each case the value of ij does not 
agree with the one that is obtained by fixing the Coulomb gauge, which 
corresponds to a = 0. Note, however, that the above gauge-independent 
results both confirm that 77 is indeed negative. 

A different point of view discussed in Ref. focus instead in the flow 
of the gauge-fixing parameter a. From the Ward identities it follows that 
the gauge-fixing parameter renormalizes as 

a = Z^ao, (97) 

which implies the /3-function 

I3 a = -lAOi. (98) 

Since at the charge fixed point we have 7a (/*)<?*) = 4 — d, the only way 
to get a fixed point to Eq. (|98(l when d £ (2,4) is to set a — 0. Due to 
the negative sign in Eq. 1)98(1 . the flow is unstable for arbitrary nonzero a. 
Thus, it is clear that stable charged fixed points can be obtained only if 
the Coulomb gauge corresponds to the fixed point of the theory in 2 < 
d < 4. For d = 4, which is the case of interest for particle physicists, 
any value of a can be chosen, since in this case 7_a(/*, g*) = 0. In Fig. 
we show a schematic flow diagram in the / — a-plane for the case where 
d E (2,4). Based on this scenario, we are led to conclude that the above 
gauge-independent results for r\ do not correspond to the infrared stable 
fixed point. This means that the correlation function l(93|l is not a good 
choice of gauge-invariant correlation function giving a gauge-independent 
value of the r\ exponent. In fact, the correlation function (f§3)> even fails to 
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u f 

Fig. 5. Schematic flow diagram in the / — a-planc. 



lead to long-range order: It can be rigorously shown that in the dimensions 
of interest the correlation function (|93H always decay to zeroP^ A better 
gauge-invariant correlation function i 



Q(r - r') = ( ip(r) exp 



le J dV'A(r") • b(r") 



^t( r ') 



where 



with 



bfr" 



W(r" 



VV(y" 



-VV(r) = S d (r) 



(99) 



(100) 



(101) 



Since the above correlation function is gauge- invariant, any gauge can be 
fixed to calculate it. The end result will be always gauge-independent. It is 
easy to see that in the Coulomb gauge 



Q{Y - r' ) | Coulomb. = (V , ( r )V ,t ( r '))|o=0- 



(102) 



Therefore in such a scenario the rj exponent corresponds precisely to the one 
we have calculated previously in the Coulomb gauge. Furthermore, it ca n 
be shown that the correlation function H99fl exhibits long-range orderp^l in 
contrast to the one given in Eq. iffify . 
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5. Renormalization group calculation at fixed dimension 
and below T c 

In the e-expansion RG approach the same RG functions are obtained no 
matter the calculation is carried out below or above T c . As dictated by the 
Ward identities the singular behavior is exactly the same above and below 
T c However, if the calculation is done in fixed dimension d = 3 and below 
T c the situation is different and the RG functions depend explicitly on the 
Ginzburg parameter re. We shall not give the details of this approach here. 
Instead, we shall concentrate on the physical aspects of this new approach 
which allows to obtain a charge fixed point at one-loop order. The interested 
reader is referred to Ref. ^|for the technical details. 

The approach we are going to discuss is not really perturbative. Actually, 
only the powers of / are being effectively counted and the powers of g are 
counted only partially. Thus, by one-loop we mean first-order in /. The 
point is that re arises in the calculations in two different ways: as the ratio 
between the masses re = m/mA and as the ratio between coupling constants, 
re 2 = g/2f. The coupling g, when it appears, is eliminated in favor of re and 
the RG flow is in this way parametrized in terms of re and /. This way of 
working is of course more complicated than more usual RG approaches but 
it has the advantage of being physically more appealing due to the explicit 
presence of re in the RG functions. In the classical Abrikosov solution of 
the GL model in an external magnetic field re appears explicitly and the 
existence of two types of superconductivity is made evident.""^ For instance, 
the slope of the magnetization curve near H C 2 is given by 

dM 1 

~dH ~ 4 7 r / 3 j4 (2re 2 - 1)' ( ' 

where Pa is the Abrikosov parameter. The above expression is singular at 
re = l/v2, which corresponds to the point separating type I from type II 
superconductivity. Such a singular behavior at k = l/v2 should be also 
visible in the GL model with a thermally fluctuating vector potential. The 
new approach introduced in Ref. IT81 makes this feature explicit in a RG 
context. As we shall see, this aspect of this new approach is crucial to 
obtain the charged fixed point at d = 3 and N — 2. 

The only RG function that is singular at K = l/v2 is 7^4: 



V2C(re)f 

™ = 24,(2^-1)3 ' (104) 
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where 

C(k) = 4k 6 + 10k 4 - 24\/2k 3 + 27k 2 + 4\/2k - 1/2. (105) 
The /3-function for k 2 is given by 

= (2 77r - 1A - ^)k 2 , (106) 

where 



K f 2k 2 + x /2k-8 
127T V2k + 1 



^ V 2 V2k, 
As before, the charged fixed point at d = 3 is determined by the condition 
7a(/*j k *) = 1- This leads to the fixed points 



/„ pb 0.3, /s» « 1.17/V2. (109) 

Note that k* is slightly above the value l/\/2 and therefore the charged 
fixed point occurs in the type II regime. 

The reason why a charged fixed point is obtained in the above analysis 
is similar to the reason why the l/A^-expansion leads to a charged fixed 
point already at order 1/N: The fixed point coupling /* is small enough 
such that a / 2 -term is strongly suppressed in the other RG functions. It is 
a large / 2 -term in j3 g that spoils at N = 2 the charged fixed point in the 
HLM theory. In order to explain why this new method is so successful, let 
us define an effective coupling / by 



7a(/» = 1- (HO) 

The above equation defines a critical line in the sense that the /3/ vanishes 
on this line. Note, however, that /3 K 2 does not vanish in general. From Eq. 
(Ill Op t we obtain 

m = 247r y~ 1)3 . (in) 
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From the above equation we see that f(n) becomes very small when n 
approaches 1 / \/2 from the right. Precisely at n = 1/ \[2 we have / = 0. This 
behavior suggests that the best approximation scheme should be one where 
the small parameter is given by An = k — 1/ \/2. Now it is easy to see that 
the e-expansion based RG fails because it effectively expands around k = 
and therefore it corresponds to the deep type I regime where the transition 
is clearly first-order. Furthermore, C(k) vanishes at K = O.QQ&/^/2 and / 
becomes very large as this value of k is approached from the left. Thus, a 
perturbation expansion around k = breaks down at n — 0.096/-\/2- There 
is an "infinite barrier" separating the deep type I from the type II regime. 
In the interval 0.096/A/2 < k < 1/^/2 the effective coupling / is negative 
and thus unphysical. The coupling / can be really small only for k > 1/ V2, 
i.e., in the type II regime. 

Note that our one-loop approximation gives only one charged fixed 
point. The tricritical fixed point is absent in this approximation. This be- 
havior also occurs in the 1 /iV-expansion where only one charged fixed point 
is found. A higher order calculation is necessary to obtain the tricritical 
fixed point. At two loops the singular behavior in n is expected to change. 
Thus, instead finding a singularity at k — l/\/2, which is the same as in the 
mean-field solution, we expect to find a singular behavior at n t ~ 0.8/-\/2, 
in agreement with Refs. Illlandll3l 

6. Concluding remarks 

In this paper we have reviewed several modern field-theoretic approaches 
in the superconducting phase transition. We have emphasized some special 
topics which are not extensively discussed in the literature. In particular, 
the scaling behavior of the continuum dual model was analysed in more 
detail than in the original publications. The duality scenario is physically 
and conceptually very important, but its scaling behavior is not yet fully 
understood. Another topic that deserves further attentio n is t he recently 
conjectured Lifshitz point-like behavior in the GL model.ESIIl Such a sce- 
nario provides an interesting possibility to understand physically the nega- 
tive value of the critical exponent r\. 
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